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Abstract 

In this work we develop on the recently suggested concept of superstatistics [C. Beck and E.G.D. 
Cohen, Physica A 322, 267 (2003)], face the problem of devising a viable way for estimating the 
correct statistics for a system in absence of sufficient knowledge of its microscopical dynamics, 
and suggest to solve it through the Maximum Entropy Principle. As an example, we deduce the 
Probability Distribution Function for velocity fluctuations in turbulent fluids, which is slightly 
different from the form suggested in [C. Beck, Phys. Rev. Lett. 87, 180601 (2001)]. 
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Nonextensive statistical mechanics (NESM) has been raising considerable interest along 
these last 15 years. However, while there is an universal consensus about its most famous 
prediction, i.e. the ubiquitous existence in nature of power-law Probability Distribution 
Functions (PDFs), the soundness of the theoretical foundation on which is based-the gener- 
alized definition of entropy [lj-has been questioned by several authors. We quote here, e.g., 
the paper j^J which dealt in particular with contradictions between theory's predictions and 
thermodynamical constraints. 

The paper jjj started from these critiques but focussed more on the information-theoretic 
aspect, suggesting that Tsallis' entropy should be regarded merely as a practical tool for 
doing predictions in presence of a reduced amount of information about the system. Indeed, 
the need of resorting to modified definitions of the informational entropy when the knowl- 
edge about the states of the system is insufficient, is well known in statistics and has been 
extensively pointed out in 141 and references therein. 

n n 

Quite recently, Sattin and Salasnich [5j]-starting from an earlier work by Beck p-and, in a 
more formal and rigorous way, Beck and Cohen ?J demonstrated, without direct reference 
to any definition of entropy, that Tsallis' statistics is just a particular case of an infinite 
class of statistics (hence the name "superstatistics"). All the elements of the class are char- 
acterized by one or more parameters, and reduce to ordinary Maxwell-Boltzmann statistics 
for particular values of the parameters. Beck and Cohen started from a model of dynamical 
system, a Brownian particle moving according to a Langevin equation, where the noise (a) 
and friction (7) terms are allowed to fluctuate, to arrive to the famous equation 

P(E) = K e-? E f((3)d(3 (1) 

J 

Here, P(E) is the PDF for the system of being in the state of "energy" E, [3 is a fluctuating 
parameter (generalized inverse temperature), which in the present formulation is function 
of a, 7, and f(/3) is the PDF for the realization of the particular value (3; e~ l3E is the usual 
Boltzmann factor. Eq. (1) is more general than suggested by the Brownian particle model. 
Indeed, it can be written for any system £ interacting with a fluctuating environment B, 
hence P(E) is a weighted average of the standard Boltzmann statistics over different re- 
alizations of the interaction between the system and its environment, quantified by f(/3). 
Notice that the existence of finite fluctuations in the environment is strictly correlated with 
finite-size-effects of the environment itself; indeed, the result that Tsallis statistics can arise 
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within the framework of a system interacting with a finite thermal bath was recently re- 
ported by Aringazin and Mazhitov 8|, and it had already been suggested much earlier by 
Plastino and Plastino 

The physical content of the theory, thus, shifts from the entropic index q to the PDF /(/?). 
Of course, the explicit expression for f(/3) must vary for any single problem, constrained 
just by some rather intuitive criteria (normalizability, etc . . . ); Beck and Cohen give several 
possible examples of functions which are potential candidate for f(/3). However, a simple 
criterion able to guide the user towards a plausible functional form for /, lacking a more de- 
tailed knowledge of the underlying microscopical details, would be very satisfying. But such 
a criterion is readily available: it is the well known Jaynes' Maximum Entropy Principle. In 
this case, the system about which we have not a proper knowledge is no longer £ but the 
environment B itself, and f{(3) is a measure of the probability of B of occupying a state in 
an abstract one- dimensional space parameterized by (5. 

Lacking any further information, the most probable realization of f{(3) will be the one that 
maximizes (Shannon) entropy S(f) = — J f In fd/3 with suitable constraints. 
In this paper we present a straightforward application of this principle to an important 
example, namely the PDF of velocity fluctuations in turbulent flows This accurate ex- 
perimental measurement is thought to be one of the strongest evidences in support of Tsallis' 
theory, since the empirical PDF appear very well matched by power-laws. In this work we 
suggest instead that the true curve can be very close numerically, but rather different in 
its analytical expression, from a power-law. Indeed, further experimental investigations of 
turbulent flows are now suggesting the existence of small deviation from pure power-laws 
experiments of Jung and Swinney, cited in |3] ) - Recent work by Aringanzin and Mazhitov 
lp| deal with the attempt of theoretically recovering the new accurate experimental distri- 
butions for fluid particle accelerations. In this work, instead, we will consider, only fluid 
velocity differences. 

As a first test case, let us suppose of knowing anything about / but the average value of the 
inverse temperature 

</?>=/ Pf(PW • (2) 
The most probable / extremizes the functional 

F = -J f\nfd{3-\J ff3d/3 . (3) 
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and the solution is the Boltzmann-like function 

/(/3) = / exp(-A/3) (4) 

(A is the Lagrange multiplier). Such an expression was obtained, by example, in from 
which the original idea of entropy maximization was taken. It reads slightly different from 
the Gamma (or \ 2 ) distribution originally devised by Beck J(| to recover power-law PDF. If 
we replace / of (jU into (JTJ) we get 

which is not a physically acceptable solution, since P(E) must be normalizable. 
In order to make a step beyond this rough scheme it is necessary to input more information 
within the model. We do it starting from the same premises as Beck's |f| but diverging at 
just the next step: infact, he defines (3 oc e r r, where e r is the energy dissipation rate of the 
fluid on microscopical scale, and r the typical energy transfer time. We note that e r r has 
units of [energy], not of [energy] -1 , and a more intuitive way of writing this relation should 
be (3 ~ (e r r) -1 . However, the inverse of a sum of squares of random variables does not yield 
a ^-distributed random variable, hence Beck is forced to re-establish the correct dimensions 
by multiplying by the constant A, with the units of [speed] 4 . Thus, a characteristic speed 
A 1 / 4 has entered the calculations, whose physical interpretation remains obscure. Of course, 
Beck was forced to do this assumption because of the sought agreement with Tsallis' theory. 
Now, we are free from this constraint, and can allow for more natural choices, although we 
must agree that a certain arbitrariness in the choice of the physically meaningful variables 
is unavoidable within this context. 

To start with, we shift from j3 parameter to T — 1//3. T, which has units of [energy] is 
a more convenient variable, as explained above. We can write the equivalent of Eq. (1) 
in terms of the new parameter; furthermore, in order to adhere to existing literature, it is 
convenient using a generalized velocity instead of energy: E = u 2 /2. Therefore Eq. (1) 
becomes 

P(«) = I dT \/S ex P 9(T) ■ (6) 



. 2vrT r \ 2T / 

Again, we follow Beck's recipe: the parameter T is written in terms of fluctuating Kolgo- 
morov velocities Ui 

° j=l,2,3 



where Tn is a constant. We stress that the r.h.s. of the previous equation and that of Eq. 
(28) in p are the same, although the l.h.s.'s are each the inverse of the other. This is due 
to our discarding of the artificial constant A. T , conversely, has a straightforward physical 
interpretation as average thermal energy. 

The PDF g(T) can be rewritten as a function of it's: g'(ui, u 2 , u 3 ) = g'(u) <-> g(T). The 
Maximum-Entropy principle imposes of extremizing 

F = -g'(u) ln(</(u)) -^( u 2 1 + u 2 2 + ul)g'(u) (8) 



u 



M 



where l/u 2 M is the Lagrange multiplier corresponding to the fact that we fix the average 
kinetic energy. The solution is, obviously, g'(u) = Cexp(— (u\ + u 2 . + m|)/mq). From here, 
reversing to T variable, 

g(T)=KT^e W (-^-) . (9) 



2T , 

The term T 1 / 2 comes from the volume element, C, K are normalization constants and we 
have rescaled velocities such that u 2 M = 2/3. Notice that Eq. (0) is a Gamma distribution 
and could be obtained straightforwardly from Eq. (JJJ by assuming from the start that the 
Ui were normal random variables, just as done by Beck. 
Let us now replace Eq. (jHJ) into (jBJ): we get 

p(u) = K' Texpt-^-^dT . (10) 



to 1 \ 2T 2T 
The explicitly normalized solution reads 



U1T U \U 




p(u) PW = --K 1 (^) \u = \l^\ (11) 



and K\ is the Bessel K function of order one. This result appears rather different from usual 
power-laws. Indeed, we chose to plot in Fig. (1) this curve together with the best fitting 
curve found by Beck for the velocity PDF: 

K«)Beck = ^r~, Z : U(a-l) (? ~ 1A > 01 ~ - 9 )- ( 12 ) 

Z -? (l + (q - l)f3C\u\ 2a ) 

(See dashed line in Fig. 1 of [qj). On the whole, the two curves match rather closely. Some 
differences appear at low u's. 

Since we have not available the original experimental data p ex pt, we cannot directly com- 
pare on them the goodness of our fit. A good insight comes, however, by plotting the 
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FIG. 1: Solid line, p(u)Beck from Eq. (|12|): dashed line, p(u)pw from present work (Eq. The 
adjustable parameters have been chosen so that the curves have the same variance. 

relative difference A = (peeck — PPw)/PBeck (Fig. 2, p PW is our solution, given in Eq. ITT|) . 
The quantity Ai = (pBeck — Pcx P t)/PBeck appears plotted in Fig. 2 of If A and Ai agree, 
Ppw — Pexpt- It appears that our fit slightly overestimates experiment at u ~ 0, but for the 
same amount as PBeck does underestimate it. On the whole, the agreement with experiment 
is remarkable. We remark that, as only adjustable parameter, we used the hypothesis that 
T is of the form (J7J) with the index ranging from one to three. The equivalent of a parameter 
used by Beck does not enter our calculations. Thus, we have realized an economy in our 
way of modelling the data. 

We call the reader's attention to the fact that, in his Fig. 1, Beck studied at once (i) 
velocity spatial differences, and (ii) accelerations. We, instead, considered only the former 
quantity, and it is possible to see that the latter cannot be reproduced by the present 
treatment, even allowing for a varying number of Kolmogorov velocities in Eq. (JJJ). This 
must be traced back to the fact that now g(T) is no longer a good weight function, but we 
must define an equivalent of T, defining the "average acceleration" of the system. 
Finally, we point to two important issues: I) we have the asymptotic trend xK\(x) — > 



x 



1/2 



exp(— x) 



x 



oo), that is, we have not power-law decay. Notwithstanding this, our 



curve nicely fits data that have been previously considered as stemming from a power- 
law PDF. The reason lies in the finite w-range sampled and therefore calls for an inherent 
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FIG. 2: The relative difference of the two curves above: A = (pBcck — PP\v)/PBeck versus u, where 
PBeck is from Eq. ()12j0 and ppw is from Eq. fjllj) . This figure should be compared with Fig. 2 of 



ambiguity in this kind of studies, related to finite experimental scans: unless one is sure of 
being investigating the true asymptotic region, one can never be completely confident about 
the fitting curve used. 

II) As reported above, all superstatistics must collapse to the single Boltzmann statistics 
when q — > 1, that is, differences between different models are at least order 0(q — 1 ). Since, 
in this case, we are dealing with a parameter not far from unity, q ~ 1.1, it is to be expected 
that any two reasonable models would give close results. (The author wishes to thank C. 
Beck for calling these two points to his attention). 

In conclusion, the use of the concept of superstatistics together with Maximum-Entropy 
principle appears to be an efficient way of estimating statistical properties in general systems 
using a minimal amount of information. 
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